Power-law behavior in the quantum-resonant evolution of the 5-kicked accelerator 
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We consider the atom-optical 5-kicked accelerator when the initial momentum distribution is symmetric. 
We demonstrate the existence of quantum-resonant dynamics, and derive analytic expressions for the system 
evolution. In particular, we consider the dynamical evolution of the momentum moments and find that all 
even-ordered momentum moments exhibit a power law growth. In the ultracold (zero-temperature) limit the 
exponent is determined by the order of the moment, whereas for a broad, thermal initial momentum distribution 
the exponent is reduced by one. To demonstrate the power law behavior explicitly we consider the evolutions 
of the second- and fourth-order momentum moments, and cumulants, for an initially Gaussian momentum 
distribution corresponding to the Maxwell-Boltzmann distribution of an ideal gas at thermal equilibrium. 
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I. INTRODUCTION 

Quantum resonances and antiresonances [2] [3] [4] [5] |6] 

EziisiEisaiiiiiEiisiiiiiiaiisiiniiisiiisiEaare 

striking signatures of ordered quantum manifestations in the 
5-kicked rotor, a paradigm system in the study of classical 
BTl and quantum chaotic dynamics E2l l23l l24l . The atom- 
optical realization of the 5-kicked rotor (more accurately de- 
noted the 5-kicked particle) has in recent years proved to 
be a fertile testing ground for quantum-chaotic phenomena 
ITT41 IT51 ri6l 1TT1 1T81 IT9l l25l 1^1 1271 1281 1291 l^Ol I3TI 1321 
[33] EH [35] [36] [371 ESI E21 |40l |H] |42] H3). An important 
variant is the 5-kicked accelerator, where the kicking poten- 
tial is aligned parallel to the local gravitational acceleration 
EU [45] [47] |48] [49] [50] EU [52] [53] [54] [55l. This gives 
rise to phenomena closely related to quantum resonances, for 
example quantum accelerator modes, and fractional quantum 
resonances |56|. Such atom-optical realizations involve sub- 
jecting a cold, dilute atomic gas to a periodically-pulsed laser 
standing wave. The amplitude, phase, periodicity, and dura- 
tion of the pulses can be controlled to a high degree of preci- 
sion, allowing theoretical predictions to be thoroughly inves- 
tigated. 

In this paper we investigate the quantum-resonant behav- 
ior e in ej is] e] noj e] [fa [i3] ni] nsi nsi uti hsj na of 

laser-driven clouds of freely-falling, laser-cooled atoms. Such 
a system closely models the 5-kicked accelerator, of which 
the (5-kicked particle is a particular case. Experimentally, the 
system dynamics are typically interrogated by measuring the 
atomic centre-of mass momentum distribution, using a time- 
of-flight technique ||46| . Although there are practical reso- 
lution limits, it is then in principle possible to determine all 
momentum moments from the measured distribution. A lim- 
ited number of moments is frequently sufficient to character- 
ize the system: the second-order momentum moment provides 
a clear signature to distinguish between quantum-resonant dy- 
namics and dynamical localization in the 5-kicked particle 
l57l . whereas the fourth-order moment is necessary to distin- 



guish between different fractional quantum resonances in the 
5-kicked accelerator at finite temperature l56l . The momen- 
tum moments therefore provide conceptually simple measures 
to characterize the effect of quantum resonant dynamics on 
the atomic centre-of-mass momentum distribution. The main 
focus of this paper is to provide fully analytical predictions 
of the momentum-moment dynamics of the atom-optical 6- 
kicked accelerator under certain physically motivated condi- 
tions. We present detailed derivations of the key analytical 
results utilized in Saunders et al. [56, 57] . 

Assuming only that the initial momentum distribution is 
symmetric, we find that the time evolution of all even mo- 
mentum moments follow power-law behaviors, and examine 
how such behavior differs between opposing limits of the ini- 
tial distribution being either extremely narrow or broad. As- 
suming a Maxwell-Boltzmann distribution for an ideal gas, 
we determine explicit analytic expressions for the time evo- 
lution of the second- and fourth-order momentum moments 
when the gas is initially in an ultracold limit, or at finite tem- 
perature. We also outline a procedure for calculating higher- 
order momentum moments for a given (but arbitrary) initial 
momentum distribution. Although we neglect two-body in- 
teractions, we note that our model can also be used to de- 
scribe Bose-Einstein condensates where the scattering length 
has been made negligibly small by tuning around a Feshbach 
resonance ll58ll59ll6CTl6Tl. 

The paper is organized as follows. In Sec. |ll] we sketch 
our theoretical model for the atom-optical 5-kicked acceler- 
ator, derive the time evolution operator from which the time 
evolution of the momentum moments can be determined, and 
derive general conditions for which antiresonance and reso- 



nance occur. In Sec. Ill we consider only momentum distribu- 
tions that are initially symmetric, and discuss the power-law 
behavior of the evolution of the even-ordered momentum mo- 
ments in the two limits of narrow and broad initial momentum 
distributions. In Sec. [TV] we consider the time evolution of the 
second- and fourth-order momentum moments and cumulants 
for momentum distributions that are initially Gaussian, that 



2 



is, corresponding to a finite temperature ideal gas. Section [Vj 
consists of the conclusions, which are then followed by five 
technical appendices. 



II. ATOM-OPTICAL .5-KICKED ACCELERATOR 

A. System Hamiltonian 

We consider a cloud of trapped and laser-cooled alkali 
atoms. The cloud is released from the trap and addressed 
by an appropriate configuration of off-resonant laser beams 
forming a vertically aligned, pulsed, sinusoidal potential. We 
neglect atomic interactions and treat the cloud as an ensemble 
of single-particle systems with Hamiltonian 



H 



2M 



+ Maz ■ 



8A 



cos(Kz)f(t). 



(1) 



Here p and z are the centre-of-mass momentum and posi- 
tion operators, Q« is the Rabi frequency, K/2 is the laser 
wavenumber along the z direction, M is the mass, A is the 
detuning and /(f) describes the periodic pulses. The parame- 
ter a is the relative acceleration between the atomic cloud and 
the optical standing-wave potential and is in general given by 
a = g — a<p, where g is the local gravitational acceleration and 
a$ can be chosen with high precision by tuning the relative 
phase of the laser beams creating the optical field [56|. When 
a is set to zero, we recover the atom-optical <5-kicked rotor 

Eil. 

We assume the pulse duration to be sufficiently short for 
the atoms to be in the Raman-Nath regime, i.e., the displace- 
ment of the atoms during the pulse is significantly less than 
the wavelength of the standing wave. Hence, to a good ap- 
proximation, the pulses may be modeled as a train of 6 func- 
tions 1 57 1 giving the quantum 5-kicked accelerator Hamilto- 
nian [!46l 



Maz - h<p d cos(Kz) V S(t - nT), (2) 

«=o 



P 

Hfikn — + 

2M 



where <pa = Q„r p /8A, t p is the pulse duration, and T is the 
pulse periodicity. For convenience we apply the unitary oper- 
ator U = exp(iMaz/n) Il52"ll53ll54l , to transform Hamiltonian 
(|2]i to the spatially periodic form: 

(p - Mat) 1 -A 
H S te = — - Hd cos(Kz) \ 5(t - nT), (3) 



2M 



n=0 



for which we can invoke Bloch theory Il57ll62l . 



B. Time evolution 



1. Transformed Floquet operator 



a Floquet (kick-to-kick time-evolution) operator. The trans- 
formed Floquet operator corresponding to the transformed 
Hamiltonian ^ is exp(-iMa 2 T 3 [3n 2 -3n + l]/6h)F„ 
where 



_ -i [hK 2 (ic+P) 2 T/2M-Ka(k+P)(2n~l)T 2 /2] M d cos(S) 



(4) 



i.e., we absorb a global phase into the definition of F„. 1 We 
have separated the momentum and position operators into dis- 
crete and continuous components l54ll : 



z = K'\2nl + 9), p = ftK(k+/3), 



(5) 



where the eigenvalues of / and k are integers, and the eigen- 
values of 6 and ft are 6 e [-n, n) and yS e [-1/2, 1 /2). 

Because the Hamiltonian Q is periodic in space, the quasi- 
momentum /3 is conserved 152H53L i.e., /? commutes with 
ll54l and only momentum eigenstates with eigenvalues differ- 
ing by integer multiples of hK are coupled [46]. Therefore, 
within a particular quasimomentum subspace, the dynamics 
are determined by 

p n (j^ _ e -i[hK 2 (k+P) 2 T/2M-Ka(k+mn-l)T 2 /2] e ifacos(e)^ ( 6 ) 

As is the case for the (5-kicked rotor or particle, quantum 
resonances or antiresonances occur when T = (Tj/2 l57l . 
where T T = 27TM/HK 2 is the Talbot time P31 (so-named in 
analogy with the Talbot effect in optics ll63l ). In this case, Eq. 
(|6]i becomes 



l-Cl(2n-l)] e -m{(t-Sl)+2tp\k gifa cos(ff) 



(7) 



where we have used exp(-i{nk 2 ) = exp(-i{nk), and intro- 
duced the dimensionless effective gravitational acceleration 
Q = KaT 2 l2n l50l . 



2. Time evolution of the momentum moments 

The time evolution of a general momentum eigenstate \k+/3) 
can be derived by consecutively applying the transformed Flo- 
quet operato r j7), the details of which are given in Appendix 
0[See Eqs. ( |A15| > and ( |A19| >]. We deduce that 



\W(t = nT)) =F„(fi)F„- 1 (fi) . . . F \(j3)\k + P) 



= ]T J hk ((jjy { i- k)x e-' mnak e-' mlm 



(8) 



./=- 



-inn0 2 ( 



\j+P>, 



where u> andx are defined in Eqs. ( A16 1 and ( A17 1, and 2 

a = (l+20){. (9) 



The 5-kicked accelerator Hamiltonian (|2]i is periodic in 
time, and the system evolution can be described in terms of 



1 The global phase is corrected slightly from that given previously [46 ] 

2 The parameter a is related to T |56 57'] by T = an/2. 
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Defining the dimensionless momentum P = k + ft = p/hK, 
it follows straightforwardly from Eq. <|8j that the m-th order 
momentum moment of the evolved state is 



(10) 



where 77 = \u)\l<pd = |v|, with 



J=0 



(ID 



We note that Eq. (jTTJ has the form of a Gauss sum, and that 
Gauss sums have applications in number theory and various 
areas of theoretical physics I64ll65ll66ll67ll . 

By extension, the evolution of the moments for an initial 
statistical mixture of momentum eigenstates, with distribution 
D k (f3), is described by 15511571 

Xl/2 00 
d P X Jl k (<i>J>l)U+P) m D k (J3). (12) 
1/2 j£L 



C. Existence of resonances and antiresonances 

1. Manifestation of resonances in n 

Quantum resonances are characterized by unbounded 
growth in the system energy, and we see that this is due to 
the constructive inteference of oscillatory terms in Eq. 



11 



They occur in the (5-kicked rotor or particle when ^77 in Eq. 
([H} can be replaced by <pdti 11571 . To achieve such resonant 
growth in 77, the oscillatory terms contained within the sum in 
Eq. (Hi must add in phase. It follows that the summand of 
v must be periodic with respect to the counting index j. The 
period Q (if present) is determined from 



>[<K/+fi)-nO' 2 +2j'G+fi 2 )] - Jniaj-af) 



implying that 



aQ-2 jSlQ - £IQ = 2A, 



(13) 



(14) 



must be fulfilled, where A is an arbitrary integer. 

The period Q should be independent of the counting index 
j, and so Q.Q must be integer. Hence, Q. is rational, and if we 
set Cl = r/s (r and s are integers with no common factors), it 
follows that the smallest possible value for Q is s. With this 



information we deduce a simpler condition from Eq. ( 13 1: 



s(a-r) = 2A. 
Eliminating a using Eq. ([9]) gives 



2A + s(r - €) 

2£s : 



(15) 



(16) 



indicating the quasimomentum values, subject to /3 e 
(-1/2, 1/2], for which quantum resonances occur. 



0.5 
-0.5 L^ 6 



oa. 



(b) 



oa. 



(c) 



oa. 



(d) 



I 

I 



1 

0.1 0.2 0.3 0.4 0.5 0.6 0.7 o.s 0.9 1 



n 



FIG. 1 : (Color online) Parameter 77 as a function of Cl and /? for 1 = 1 
at (a) 15 kicks (b) 30 kicks, and (c) 60 kicks. An antiresonance of 
order 2N S , with Cl = r/s, will initially mimic exactly resonant growth 
for N s s kicks. The diminishing finite size of the resonances in (a) - 
(c) is due to nearby high-order antiresonances which become appar- 
ent as n increases. Panel (d) shows the location of those resonances 
(see Eq. l |16[ l) for Cl = l/s, where ^ = 1, 2, . . . , 20. The resonances 
become infinite in number as Cl — > 0, but for Cl = there are only 
three resonances. 



2. Manifestation of antiresonances in 77 

An antiresonance, in the strictest sense, is when the state 
periodically cycles back to its initial condition [57|. This 
can only occur when ^77 is periodic in «, implying that 
exp(/7r[a7 - Q/' 2 ]) should be oscillatory, i.e., 



piTTlaU+QyCXf+ljQ+Q 2 )] _ _ Maj-af) _ in(aj-Cij 2 + l) 



(17) 



Proceeding along similar lines as for the resonant case (see 



Sec. II C 1 1 this implies that Q — s, and that 
s{a-r) =2A + 1, 



(18) 



where A is an arbitrary integer. Condition (I81 ensures that 
77 = when n is an even multiple of s. 

For a given value of Q = rj s, higher-order antiresonances 
also exist for different values of/? [57|. These are described 
by the condition 



e 



i"la(J+Q)-(r/s)(j 2 +2jQ+Q 2 )] _JnN r /N s i n [ a j-(r/.s)j 2 ] 



(19) 



^ e i7Tla(j+N s Q)-(r/ s )(j 2 +2jN s Q+N 2 Q 2 )] _ _ g 6r[aj-(r/s)/] _ 

As above, it follows that Q = s; hence, the revival period of 



tion of Eq. ( 18 1, i.e., 



the initial state is 2N s s. Condition ( 19 1 implies a generaliza 



s(a - r) = 2A + 



N, 



(20) 
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Antiresonances with revival periods 2N s s therefore occur for 
quasimomentum values (subject to ft e (-1/2, 1 /2]) 

2A + N r /N s + s(r - €) 



lis 



(21) 



How both resonances and antiresonances are manifest in rj, 
and their dependence on Q and ft, is illustrated in Fig.[T] This 
demonstrates that the resonances are dense in the parameter 
space, but diminish in strength for increasing s. Furthermore, 
resonances become more sharply defined for increasing n as 
nearby high-order antiresonances become more important. 



III. TIME EVOLUTION OF THE MOMENTUM 
MOMENTS 

A. General symmetric momentum distributions 

1. Physical motivation 

In the atom-optical kicked systems under consideration, the 
atomic gas is initially held in a harmonic trapping potential. 
The ground state of the system is therefore spatially symmet- 
ric, both in the presence and absence of significant interatomic 
interactions. We consider the momentum moment dynam- 
ics for the case of symmetric initial momentum distributions 
D(P), such that D(P) = D(—P), In particular, we provide a 
detailed analytic investigation of the second- and fourth-order 
momentum moments. These moments are most relevant as 
signatures of dynamics in the atom optical delta-kicked par- 
ticle 11571 and accelerator |56|. The second-order momentum 
moment [m — 2 in Eq. ( 12 1] is proportional to the mean ki- 
netic energy, and the fourth-order moment contains informa- 



tion about the degree to which the distribution is peaked. 



2. Consequences of symmetry 

We consider the initial state to be an incoherent ensemble 
of momentum eigenstates, i.e., the momentum representation 
of the density operator is diagonal. For a symmetric initial 
momentum distribution, the initial population of state | -k—ft) 
is always equal to the initial population of state \k + ft). Using 
Eq. (JTOj, and explicitly stating the /^-dependence of rj(ft), for 
an initial state | — k — ft) the momentum moments evolve as 



</*">„ = (-!)" 



2 jj_ k (cf> d ri(-PW + Pf 



(22) 



]=- 



where we have relabelled the summation (setting j to —j) and 
used that /_/jc) = (- 1 ) j Jj(x) J68] • 

Comparing Eqs. (22 1 and ( 10 1, we see that when J](-/3) - 
Tj(fi), it follows that <P m+1 >„ for an initial state | - k - ft) is 
equal to -(P 2m+1 )„ for an initial state \k + ft). The momentum 
moment for a statistical mixture is the normalized sum of the 
moments of the individual component states. Consequently, 
for initially symmetric momentum distributions, the odd mo- 
ments are identically zero whenever t](ft) is in general equal 
to t](-ft), i.e., for all values of the quasimomentum ft. 



3. Invariance ofrj 

The invariance of rj(J3) upon changing ft to -ft is a useful 
property which is often satisfied. In particular, for rational 



Q. = r/s it is valid at resonance [see Eq. (15 1] when n is a 



multiple of s, and at antiresonance when s is an even multiple 
of s. Here we derive these results. 



To begin we write Eq. (Hi, changing the summation index 
j to n - j, i.e., 



V n+l (p) = V e im+Wn-j)-(r/s)(n--2nj+j 2 )] 
n 

^ e iii(\+2l3)tn-( r l s )n 2 e ij<[-(l+2fi)tj+2(.r/s)nj-(r/s)f]^ 
J=0 



(23) 



where now we use Eq. ( 23 i to describe the n + 1th kick. Noting 
that e ml i = e~ ml i and taking n to be a multiple of s (i.e., n = 
ts) yields 

TS 

v TS+x (P) = e'V+W"-^ e im-2PW-(r/s)j 2 ]^ (24) 
;=o 

Hence, as r/ TS+1 (f3) = |v„+iQ3)|, it follows that t] TS +\(J3) = 
^ts+i ( — /?) m general, whenever £2 = r/s. If the initial momen- 
tum distribution is symmetric, all odd moments will be zero 
at these times; for integer £2, this means the odd moments will 
always be zero. 



If n = ts - 1, we deduce from Eqs. (23 1 and (24i that 



v TS (J3) =^e inl(1+2l3)eHf,s)j2] -^(i+WiT-rjT 2 



7=0 

in(l+2/3){sT-rsr 



(25) 



Hence, it follows that rj TS (J3) = rj TS (-/3) for the values of yS 
satisfying (1 + 2j3)EsT - rsT 2 = 2A (i.e., for ft an even inte- 
ger). This condition is fulfilled whenever Eq. (JT3J or Eq. ( p"8j ) 
(for even f) hold, i.e., for values of ft where resonances and 
antiresonances are supported. 

Finally, if the initial momentum distribution is symmetric, 
all odd momentum moments will be zero when n is a multiple 
of s, for resonant evolution, and when n is an even multiple 
of s, for antiresonant evolution. For this reason, we consider 
only even momentum moments in the remainder of this work. 



4. Time evolution of the even momentum moments 

The time evolution of the even momentum moments is 
given by 

Xl/2 00 
dft V jj(<p d ri)U + k+ft) 2m D k Q3), (26) 
1/2 j£L 
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where we have shifted the index j in Eq. ( 12 1 by k. Note that 
[by Eq. (JTTJ] rj(J3) - rj(k + /J). Making the change of variables 
P — k + p and binomially expanding (j + P) 2m implies 

dPD{P) V 2( ™) J J^j 2hp2(m ~ h) ' 

j=-oc h=0 ^ I 

(27) 

where D(P) = Dk(J3) and we have eliminated odd powers of 



j, because they sum over j to give zero [see Eq. (B5 1]. Using 
the general form of the summation of Bessel functions over 



even powers of j [see Eq. (B7 1], we find that 

{P 2 "\ = {P 2 "') + J dPD(P) 2h \R 2h (<f>dTl)P m ~ h \ 

(28) 

where R 2/, is a 2/ith degree even polynomial in (p^r]. 

Considering the two simplest even momentum moments 



[see Eq. (B8 1], Eq. (28 1 for m — 1 becomes 

2 



dPD(P)r, 2 , 

CO 



(29) 



and for m = 2 becomes 

&4\ _ itA 



L dpD( 4 



(30) 



B. Ultracold limit 

1. Time evolution of the moments at zero temperature 

We consider the most trivial example of a symmetric initial 
momentum distribution, i.e., we choose D(P) = 6(P) such that 
all the atoms are initially in the P = state. This describes an 
ideal zero-temperature gas, which we refer to as the ultracold 
limit. In this case ft = 0; hence a — I, and the resonance 



condition (15 1 becomes s(r - £) — 2A. Equation (28i then 
simplifies to 



(P 2m )„ = R2m(<f> d n), 



(31) 



i.e., the growth of (P 2m )„ scales to leading order as (0df/) 2m , 
as shown for integer D. in Fig. [2] In particular, we obtain 

/S2\ 2 itA\ 4,^2 /o-n 

<^ )n = yiT, (P )n = -g-77 + . (32) 
for m — 1 and m = 2, respectively. 



2. Integer values of Q. 
The case where Q = r/s is an integer is equivalent to choos- 



ing s = 1. The amplitude of Eq. (Hi is then given by Eq. 
( |A22| > as 



sin(n[^-r]7r/2) 



sinfl^ - r]n/2) 



(33) 




FIG. 2: (Color online) The first seven even momentum moments of 
a <5-kicked atomic cloud in the ultracold regime as given by Eq. |3TJ 
for integer Q, with the resonance condition l |15| l satisfied by choosing 
any appropriate value of I, We choose d> d = O.Sn as corresponding 
to an experimentally typical driving strength [46|. The growth tends 
asymptotically to a power-law behavior with the power given by the 
order of the moment. 



which identifies r\ as the modulus of Chebyshev's polynomial 
of the second kind l69ll70l . t/„_i (cos([^- r]n/2)). Imposing 
the resonance condition s(r - £) — 2A with 5 = 1 then gives 



rj - n. We therefore deduce from Eq. ( 32 1 that 



(34) 



Similarly, imposing the antiresonance condition p0| ) im- 
plies t - r — 2A + N r /N s . Hence, in the antiresonant case, 



77= C/„_ 1 ((-l) A cos(7rM./2^)), 



(35) 



and the antiresonant evolutions of the second- and fourth- 
order moments for integer Q are given by 



(P 2 ) n =<P 2 > + y t/ 2 _! (cos(nN r /2N s )), 

O id 

(P 4 ),, =<P 4 > + U li (cos(7rN r /2N s )) (36) 

+ ^t/ 2 _ I (cos(^./2^V s )). 
As an example, we consider the simplest antiresonance with 



N r = N s = 1, so that U n -i = (-1)" sin(«7r/2). From Eq. g§ 
we find that 



1 z 

(P\ ^o + y^-ir' + i], 

(P 4 )„ =^ 4 )o + (^ + f)[(-D" +1 + l]' 



(37) 



demonstrating that the second- and fourth-order momentum 
moments oscillate with a period of two kicks. 
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Rational values offl 



From the definition of rj [see Eq. (JTTJ], it follows that 



If Q. — r/s, then, for n > s, v may be divided into identi- 
cal summations with a remainder term. Hence, assuming the 



resonance condition (15 1 to be fulfilled, setting n — ts + A 
gives 



'2> 

7=0 



in[tj-(r/s)j 2 ] 



A-l 

^ e &r[fy-(r/*)j 2 ]_ 
7=0 



(38) 



It is possible to rewrite the first summation using the reci- 
procity formula lITTIl : 



ICI-l 

I 

7=0 



^ e m[(Aj 2 +BJ)/C] 



\A\-l 



m[{\AC\-B 2 )IAAC] 



2>" 



m[(Cf+Bj)/A] 



7=0 



(39) 

where A, B, and C are integers such that AC - B is even. This 
restriction is in fact exactly equivalent to the resonance condi- 



tion (15 1. Consequently, 

r- |r|-l \-\ 

y _ T Ji e iT(t 2 s/4r-r S ) ^ { B[(lM)/f] + V ^H'/'l/l, (4Q) 

' r 7=0 7=0 

We can then impose /I = because, for large values of t un- 
der conditions of quantum resonance, the second summation 



in Eq. ( 40 1 is only a small fluctuation in v. Taking the simplest 



rational value of Q = l/s, Eq. (40 » then simplifies consider- 
ably to give 



7] = T^fs. 



(41) 



Substituting Eq. (41 1 into Eq. (32i, and letting t = n/s, yields 

is 2 



yd 2 

2s 



"d 4 , Vd 2 



(42) 



where it should be noted that these expressions are strictly 
valid only when n is a multiple of s. Hence, in addition to 
the conventional quantum resonances associated with the 5- 
kicked rotor/particle |57| (i.e., when Q — 0) there are numer- 
ous fractional quantum resonances observable in the <5-kicked 
accelerator when Q takes rational values. The manifestation 
of such fractional quantum resonances in an atom-optical con- 
text at finite temperature we address in detail in a companion 
paper l56l . 



C. Finite-width momentum distributions 

1. General Form of rj 1 ' 1 



From Sec. IIIB it is apparent that elucidation of the time 



evolution of the momentum moments is dependent on knowl- 
edge of the behavior of rj. It is not feasible to evaluate rj for all 
cases. However, it is possible to consider its general behavior 
within certain approximations. 



where 



n-l 



V e btaj-af) 
7=0 



2<l 



n-l 



. V e <Mni,)+/;ca,)-nj;)-/sG<j;)] j 
W?=o 



(43) 



F(j q ) = Yffjr - nfi 

r=l 

9 



(44a) 
(44b) 



r=l 



For brevity, we have used the convenient shorthand Zj j' 
to denote the 2q sums, with j 9 = (j\,.-.,j q ) and j' = 
(J' v . . . , j' q ). Substituting Eq. d43l into Eq. d28]i then yields 



m in, \ h n-l 

"=i v ' 9=1 j,j:,=o 

Xco 
dPD{P)P 2{m - h) e m >" ) - G ^ ) \ 
oo 



(45) 



where we have expanded R2h(<f>d r l) = 2jo=i a q(<f'd r i) 2q '■> an d re_ 
placed Bby k + B = P. 



2. Large finite-width limit 



The integrands in Eq. ( |45] l are in general oscillatory and, 
by the method of stationary phase, under many circumstances 
give a negligible contribution to the momentum momen- 
tum evolution. We first separate out the integrals with non- 
oscillatory integrands. In particular, if G(j 9 ) = G(j' q ), then 
F(jq) ~ FQ' q ) = 2QpQ q ,j' q ), where p(j q ,j' q ) is an integer given 
by 3 



1 L f 



(46) 



x=l y=x+l 



3 A simple example is for q = 2 where we have the constraint j\ + jo = 
/j + f~. Considering (ji + y'2) 2 = (j'J + j' 2 ) 2 leads to the conclusion that 



j 2 i+jl-fi 2 -j' 2 2 = 2(f 1 f 2 -jij2). 
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Hence, partitioning Eq. (45 1, and using Eq. (46 1, yields 



h 

x Z fl ^ 9 z 

9=1 



Ga,)=c(j;,) 



(47) 



ZdS-rf E 



h=\ N 7 ^=1 



>[^gv-ju)] 



J<7 Jf/ 



J" 



X dPD(P)P l{m - n, e 



2(m-h) mP[GQ q )-GQ')] 



In Eq. (47 1, the number of terms in the first summation over 
j g and is equal to the number of different ways G(j 9 ) = G(j' CJ ) 
can be satisfied. This is given by a degree 2q - 1 polynomial 
in «, i.e., 



2q- 1 

S 24-1(1) = 2^ M r - 



(48) 



The coefficients b r can readily be computed for a given q, as 
outlined in Appendix [C] Note also that, for integer values of 
Q, exp(/27r£2p(j 9 , j')) collapses to unity, independent of the 
value of p(j q , j' ). Hence, for integer Q, 



E > 



(49) 



and substituting Eq. (49 1 into Eq. (47 1 yields 



+E( 2 2 I)l>? E «• 

A=l X '9=1 

00 



.fatfa,)-^;)] 



It is illustrative to rewrite 



(50) 



(51) 



from which it is apparent that Eq. ( [BT} has a maximum period 
in P of We can therefore consider initial distributions 
D(P) with characteristic width » l/£ to be "broad," meaning 



that the oscillatory terms in the integrals of Eq. ( 50 1 cause 
them to average essentially to zero. This leaves 



m 1 h 

h=\ \ Zn ' q=l 



(52) 



where the leading order is n 2m_1 . Therefore, for sufficiently 
broad distributions, we expect the 2mth-order moment to scale 
with the number of kicks n as n 2m ~ l for integer values of O. 

For Q = r/s, the terms take s different values. 

If these values to some extent add in phase, then we expect 
the growth of the 2mth-order moment to have the same power 
law as for integer values of Q., but with a constant coefficient 
that may depend on Q.. As described in a companion paper 
ll56ll . through a combination of numerical and analytical in- 
vestigation we have found that, for a broad Gaussian initial 
momentum distribution, the growth of the second-order mo- 
ment appears to be largely independent of the value of Q. The 
growth of the fourth-order moment, however, is generally cu- 
bic with n for rational £2 - r/s, but with a constant coefficient 
that diminishes with increasing s. 



D. Cumulants 

The 2mth-order moment also includes information regard- 
ing all lower-order moments, so observations of the moments 
alone do not explicitly isolate effects of different order. How- 
ever, the moments {P'") may be used to construct cumulants 
«F"» EU |72] EIC! EH, which ^ independent quantities. 
The first-, second-, third- and fourth-order cumulants are the 
mean, variance, skew, and kurtosis, respectively, and are given 
in terms of the moments by: 



<(P» =(P), 



«P 2 )) 



--(P 2 ) ■ 



(py 



((p )) =(n - XP)(p z ) + 2(py\ 

((P 4 )} =<p> 4 - 3<P 2 ) 2 + \2(P) 2 (P 2 ) - 6(P) 4 . 



(53a) 
(53b) 
(53c) 
(53d) 



The skew quantifies the asymmetry of the distribution about 
the mean, and the kurtosis quantifies the degree to which 
the distribution is peaked. For example, a Gaussian (or 6- 
function, which can be defined as a zero-variance limit of a 
Gaussian) has kurtosis = 0, whereas a distribution which is 
more sharply peaked or cusp-like has positive kurtosis, and a 
distribution which is more "blunt" has negative kurtosis. 



As discussed in Sec. Ill A we are largely considering sym- 
metric distributions, where all odd moments, and therefore 
cumulants, are zero. In this instance, the two lowest-order 
non-zero cumulants are given by: 



((p 1 )} = <p z >, «p 4 » = (py - 3(P 2 } 2 . 



(54) 



The evolutions of the second- and fourth-order moments in 
the ultracold regime [D(P) = 6(P)] for Q = r/s are described 
by Eq. ^2\. From this the lowest even cumulants, as defined 



in Eq. ( |54| >, readily follow: 



'd__2 



2s 



2 

— n ■ 
2s 



8s 2 ' 



(55) 



Hence, the leading-order power-law behaviour manifest in the 
lowest two even moments is also observed in the lowest two 
even cumulants. We have investigated numerically the be- 
haviour of all even cumulants up to m - 100, which also 
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exhibit the same leading mth-order power-law behaviour as 
the rath moments. The scaling laws therefore do not appear to 
be simply manifestations of lower-order effects. 



Although the cumulant definitions of Eq. ( 53 i are formally 



convenient, determining the moments from an experimentally 
measured momentum distribution and then combining them 
into cumulants will cause experimental errors to accumulate. 
Alternative, more direct assessment of the "width" and "point - 
edness" of the measured distribution, providing it yields es- 
sentially the same information on the system dynamics, may 
turn out to be experimentally more convenient. 



IV. FINITE TEMPERATURE IDEAL GAS 



A. Overview 



In Sec. lIHC2l we determined that for an initial momentum 
distribution that is symmetric and in some sense sufficiently 
broad, the 2rath-order momentum moments are expected to 
grow as n 2 '" -1 , where n is the kick number. That is, the power- 
law growth is in general reduced by one compared to the ul- 



tracold limit (Sec. Ill B I. We can illustrate this with analytical 
expressions for the second- and fourth-order moments and cu- 
mulants when the initial momentum distribution is that of an 
ideal gas in thermal equilibrium, i.e., 



D(P) = 



1 



V27TW 2 



exp(-P 2 /2w 2 ). 



(56) 



Here w is the standard deviation and the corresponding Boltz- 
mann temperature is given by T K = Ti 2 K 2 w 2 1 Mk B , where k B 
is Boltzmann's constant. 

Here we choose to constrain Q. to take integer values r, that 
is we do not consider the fractional quantum resonances de- 
rived for Q. = 1 Is in the ultracold limit in Sec. |IIIB 3| We 
note, however, that when considering the manifestation of 
fractional quantum resonant dynamics in a finite temperature 
ideal gas, the behaviour of the fourth moment distinguishes 
clearly between different kinds of fractional quantum reso- 
nances, whereas the behaviour of the second moment does not 
l56l . This provides an additional, specific motivation for bet- 
ter understanding the dynamics of the fourth-order momen- 
tum moment. 

When comparing with experiment, we note that our ana- 
lytical results implicitly assume that the momentum distri- 
bution can be determined with perfect precision. Obviously 
this cannot be exactly fulfilled experimentally (although Gus- 
tavsson et al. Il76l report measuring the momentum distribu- 
tion of cold caesium atoms with a resolution equivalent to 
O.OSM'). 4 Our numerical investigations indicate that, at least 
for the second- and fourth-order moments, the qualitative be- 
haviour of the moments is not very sensitive to coarsening of 
the momentum resolution. Also note that, as quantum reso- 
nant dynamics are associated with the momentum of a portion 



of the atomic cloud increasing ballistically, as time progresses 
a proportion of the atoms will no longer be in the Raman- 
Nath regime. Experimentally, we therefore expect the pre- 
dicted power-law behaviour to be followed for a finite time 
only, as dictated by details of the experimental configuration. 



B. Calculation of the Moments 



The general procedure is outlined in Appendix IdJ starting 



with Eqs. (29 1 and (30 1, we evaluate 77, expand the resulting 



expression in terms of cosines, and then integrate. Using Eqs. 



(29j) and Q to determine Eqs. p35) l and < |D9) , we find that 

d> 2 n_1 

(P\ =w 2 + fn + <f>j Yj,-D q(e - r \n - q)e- 2 « 2e2 ™\ (57) 



7=1 



and 



,A4, , 4 , §d,~ 2 , i\ , , 1 2,2 

(™) =3w H (2n + l)n H n + iw <p d n 

8 2 

5/ n-l 

+ s Jd_ Y^(_iy(e-r) e -2q 2 e 2 x 2 w 2 

^ <7=1 

2n 

q 3 - 2nq 2 -q + y(2« 2 + 1) 

1 i4 2n— 2 

ll ^(_ 1 )«(f-r) e -2 9 W W 2 



q* 2 , (12« 2 - 1) 2n 2 

Y~ 2nq + — 3 — q + Y ~ 

+ ^'f J (-ir ie - r \n-q)e- 2 ^ W 

q=\ 
n-l 

+ 6(/> 2 d Yj{-Vf (f: - r \n - q)w 2 (\ - 4w 2 q 2 e 2 x 2 ) 



(58) 



xe- 2 ^ eW . 

The timescale over which the resonant growth of the second- 
order moment (where ( - Q. is even) switches from the ul- 
tracold regime to the large finite-temperature regime can be 



deduced from Eq. (57 1. A substantially diminished growth of 
the second-order moment occurs when q 2 £ 2 n 2 w 2 » 1, so we 
take hq = l/lmv to be the transition time. Assuming that w 
is sufficiently large that «g < 1, the exponential terms in Eqs. 
p7j ) and p8j ) become vanishingly small, and we find that 



{P 2 ) n =w 2 + ^n, 



(59a) 



Note that TiK is equal to two photon recoils. 



<P 4 )„ =3w 4 + -^(2n 2 + 1)« + -f (6w 2 + 1)«. (59b) 
8 2 

The transition time characterizing the transition from quartic 
to cubic growth of the fourth-order moment is not straightfor- 
ward to define. However, from Eq. ( |58] l and Fig. [3] it is clear 
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FIG. 3: (Color online) Second-order (a) and fourth-order (b) mo- 
mentum moments yielded by the quantum-resonant evolution of an 
atomic cloud, as propagated by Eq. |7). We set fl = (i.e., no net 
acceleration), T = T T (C = 2), and choose fa = 0.8/r |46 |. The initial 
momentum distributions are Gaussian [see Eq. j56}] with standard 
deviations w given by (*) w = 1/8, (A) w = 1/32, ( ) w = 1/128, (o) 
w = 1/512. The solid lines correspond to analytic results given by 
Eqs. j57} and l |58} , and the symbols correspond to data from a Monte 
Carlo simulation, the details of which are discussed elsewhere (57). 
The dashed-dotted lines indicate the ultracold limit [see Eq. {55} with 
s = 1], and the dashed lines are a lower bound limit for a broad dis- 
tribution [see Eq. ( |59fr with w set = 0]. The vertical dotted lines in 
(a) indicate n c , = l/tnw for w = 1/8, 1/32, 1/128, i.e., the discrete 
time characterizing when growth in {P 1 } deviates from quadratic to 
linear. 



that this occurs over the same timescale as for the second- 
order moment. In previous work we considered a symmet- 
ric initial momentum distribution that is uniform over a finite 
range (less than |57l ; the same qualitative behaviour is 
then manifest as in Fig. with comparable timescales. 

C. Momentum cumulants and the thermal limit 



As discussed in Sec |III D( the momentum moments are in 
general not independent quantities and include lower-order 
correlations. It is therefore instructive to consider the second- 
and fourth-order cumulants derived from a Gaussian initial 
momentum distribution [Eq. (56i]. From Eqs. (54i and (59 1 
it follows that {{P 2 ))„ = {P 2 )„ and 



((P 4 })n 



(2« 2 -6n + \)n + 



<? d n 



(60) 



Hence, as in the ultracold limit [Eq. (55 1], the power laws de- 
rived for the momentum moments in the thermal regime also 
appear to hold for the cumulants; the long-term cubic power- 
law growth of {{P 4 ))„ is illustrated in Fig. |4] As in the ultra- 




FIG. 4: (Color online) Cube-roots of the fourth-order momentum 
cumulants of a resonantly kicked atomic cloud evolved according to 
Eq. jTJ. Parameters are SI = 0, T = T T (I = 2), and we choose 
fa = O.Sn |46|. We consider Gaussian initial momentum distri- 
butions with standard deviation (*) w = 1/8, (A) w = 1/32, (□) 
w = 1/128, (o) w = 1/512. After a certain number of kicks, the 
growth becomes linear, i.e., the fourth-order cumulant grows cubi- 
cally. The inset shows the ideal (assuming perfect 5-kicks) long-time 
behavior. 



only shown this for the second- and fourth-order cumulants, 
we expect it to be generally true. 

Note also that both the moments and the change in the 
fourth-order moment (i.e., (P 4 )„ - (P 4 )o) formally diverge for 
large w, whereas this is not the case for the changes in the cu- 
mulants from their initial values. It is therefore in the cumu- 
lants, rather than the moments, that one can speak of a well- 
defined thermal limit 1561. 



cold limit (Sec. HID i, the scaling laws are therefore not sim- 
ply manifestations of lower-order effects. Although we have 



V. CONCLUSIONS 

We have studied the quantum-resonant dynamics of the 
atom-optical quantum 5-kicked accelerator, a fundamental 
system in the study of quantum chaos which reduces to the 
<5-kicked particle (the atom-optical realization of the quantum 
<5-kicked rotor) when the (tunable) effective acceleration pa- 
rameter a = 0. In the (5-kicked particle, quantum resonant 
antiresonant dynamics result when the periodicity of repeated 
off-resonant laser pulses is a half-integer multiple of the the 
Talbot time T T = InM/hK 2 . We have found that, for such 
pulse periodicities, fractional quantum resonances occur when 
the rescaled effective acceleration parameter Q = aKT 1 /2n 
takes non-integer rational values. 

We have focussed on the dynamics of the momentum mo- 
ments of the atomic centre-of-mass momentum distribution. 
The moment evolutions are useful for characterizing the dif- 
ferent kinds of dynamic behavior in the atom-optical <5-kicked 
accelerator. We have found explicitly that, in the ultracold 
limit, quantum resonant dynamics cause all even-ordered mo- 
mentum moments to exhibit a power-law growth, with an ex- 
ponent equal to the order of the moment. For fractional quan- 
tum resonances the moment growth follows the same power 
laws, but with coefficients that reduce with increasing size of 
the denominator of LI. We have argued that, for a sufficiently 
broad initial atomic centre-of-mass momentum distribution 
undergoing quantum resonant, or fractional quantum resonant 
evolution, the exponent of the power law growth of the even- 
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ordered moments will generally be reduced by one. In every 
case considered, we have found that the long-term momen- 
tum cumulant dynamics also display power-law growth with 
the same exponent, meaning the scaling laws in the growth 
of the moments are not simply manifestations of lower-order 
effects. 

We have illustrated a reduction in the power-law exponent 
by calculating explicit expressions for the dynamics of the 
second- and fourth-order momentum moments for a finite tim- 
perature ideal gas initially at thermal equilibrium, i.e., using a 
Maxwell-Boltzmann distribution to describe the initial state. 
These are the two simplest moments to calculate, and are 
physically motivated in that they provide directly useful in- 
formation on the dynamical behavior of the atomic cloud: the 
second-order moment provides a signature to distinguish be- 
tween dynamical localization and quantum-resonant dynam- 
ics in the atom-optical <S-kicked particle [57| and accelera- 
tor, and the fourth-order momentum moment distinguishes be- 
tween different classes of fractional resonant dynamics in the 
quantum 5-kicked accelerator ll56ll . 
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APPENDIX A: MOMENTUM-STATE TIME EVOLUTION 
1. Overview 

In this appendix we determine the time evolution of the 
eigenstate \k + /3) by applying the Floquet operator ([7]) consec- 
utively for n kicks. Our derivation follows closely the equiva- 
lent calculation for the ^-kicked rotor (where Q = 0) ll57ll . 

2. Application of the transformed Floquet operators 

The effect of applying the nth transformed Floquet operator 
(j7]l to the momentum eigenstate \k + /3) is 

F„(j3)\k + P) = j dzF n (J5)\z)(z\k+P) 

=e W+M-K7n]P (Al) 

x j dze-^e^^zXzlk+fi), 

where we have defined 

Ky n = n[(l+2P){-Q.(2n-1)]. (A2) 



Using exp(-ikK~y„)\z) = exp(if3Ky n ), an d (z \k + /3> = 

i that 



■\fKj27T exp[i(k + fi)Kz], we deduce from Eq. ( Al 
l2n e J 



F n (J3)\k+P) 



dz\z + y n )e' 



i(k+P)Kz 



(A3) 



where we have used {z\k + /3> = ^jK/l^e 1 ^^ . 

The combined effect of the nth and (n + l)th transformed 
Floquet operators can be determined similarly to obtain 

F n+x (p)Fn(m+P) = ^-e i2 ^ l+ W J" dz\z + y n + r« + i> 

x e i(k+P)Kz e iMcos(Kz)+co S (Kz+Ky,,)] 

(A4) 

Hence, the full time evolution of the initial momentum eigen- 
state, governed by f„(J3) = F„(/3)F„^(J3) . . . (/3), is 



f„(B)\k + 0) =J— e '"^ (1+ « f 
V 2n 



I dz Z + ^ J n < 

J n'=l 

Xe» + ^expU d 



cos(Kz) 



(A5) 



+ f j coslKz + f j Kyj)]) 

;=i V /=i /J/ 



3. Spatial representation 

To determine the matrix element (z|!F„(/5)|£ + /3), we project 



Eq. (A5 i onto \z), i.e., 



V 2n 





( 




x exp 




COS 






V 


n-l 






+ 1 


cos 


Kz- 









Kz-^Ky* 



y n'=\ 



j'=\ ) 



Using Eq. l |A2) , we can readily evaluate the sum 
Ky„- = nn [(1 + 2/3)1 - Qn] , 



and, defining 



qj = (n-j)n[(l+2fi)£-Qin + ./)]. 



(A6) 



(A7) 



(A8) 



Eq. ( A6 1 simplifies to 



(z\T„(J3)\k +P) = a/ -— e ™^M e ^^-^ 
" 2n 



I n-l 



x exp 



i(f> d ^ cos(Kz - qj) 



\ 7=0 



(A9) 
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Finally, defining ^ = Y/j=o COS( lj an d C = Z" = d sm< 77> Eq. 
(|A9|> further simplifies to 



x g'fcf cos(A:z)+i0jf sin(Tft) 



fe Evaluation of o) and x 

It is convenient to define 5 

0</v = we'*, (A16) 



i.e., 



4. Probability amplitudes 



a. General form of the probability amplitudes 



7=0 



Substituting Eq. ( A8 i into Eq. ( A17 1, we obtain 



(A17) 



Invoking Bessel function expansions [57|,Eq. (AlOlcanbe 
recast in the form 



(Z\f n (j3)\k +0)=J L^WJOMfzr*) 

V 2n 



x J] e iiKz Jj{cS)e ijx , 



(All) 



where to andx are real and defined by we'* = <pd{i^ + 0- To 
evaluate the matrix element 

{j+f?\r n q3)\k+P)= yj^f dze- iW ' )K \z\F n (J3)\k+P), 

(A12) 



we substitute Eq. ( Al 1 1 into Eq. ( A12 1 to yield 



{j+P\T„(f!))\k+P) = e -i™<*k e -in 2 n(k+[i)Sl e -innl3il ^ j ^ 

j'=-oo 



x e - 



,K_ 
2n 



dze 



i(k+f-M-p)Kz 



(A13) 



where a = (1 +2/3)(. 

We now consider the general expansion 



mt = bD> = f„(/3)\k +P) = J] c kj (/3,nT)\j+f3), (A14) 



where the probability amplitudes Ckj(f3,nT) are given by 
q 7 (/3, nT)5(p- p) = {j+[3'\f n (J3)\k+P). Evaluating the Fourier 
integral in Eq. ( |A13| l, it follows that 



c kj (f3,nT) = J j _ k (co)e iU - k)x e- imak e- in2n( - k+m e- i "^ 2e , (A15) 



where the normalisation condition £ J =oa \ckj(fi, nT)\ 2 = 1 is 
satisfied. Equations ( A14 1 and ( A15 1 describe the time evolu- 
tion of the initial momentum eigenstate \k + f3). However, it 
still remains to evaluate w and^- 



,,— mificn— On 2 ) ^ ^ 



in(aj-Clj 2 ) 



(A18) 



7=0 



Absorbing the global phase, we define 

n-l 

y _ V e «(«7'-n/) ) 
7=0 



(A19) 



and hence the magnitude 77 = |v[ = |v| = M/^. It is per- 
haps natural to think of a> in we'* as being positive, although 



we find it more convenient in Eq. ( A22 1 to allow oj to take 
negative values. 



c. Integer values of CI 



The Gauss sum in Eq. ( |A19[ > can be evaluated analytically 
in some particular cases. Here we illustrate this for integer 



Q = r. For this choice of Q, Eq. ( A19 1 reduces to a geometric 
sum, 



n-l 

7=0 



(A20) 



which can be evaluated to give 

1 _ ^ _ J(n _ l)(a _ r)n/2 S m(n[a-r]7T/2) 



v = 



sin([a - r\n/2) 



Referring to Eqs. ( A17 1 and ( A18 1, we can now set 

w sin(n[o' - r\n/2) 
4>d sin([a - r]n/2) ' 



(A21) 



(A22) 



and^ = [1 - (a - r)(n + l)]7r/2. Note that Eq. \A22\ can take 
negative values, and hence the absolute value must be taken to 
obtain 77, if desired. 



5 The parameter v in this paper is related to a previously defined parameter fi 
1 57 1 according to v = i/i*. 
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APPENDIX B: SUMMATIONS OF BESSEL FUNCTIONS 



We first state the recurrence relation 

In 



Jn-l(x) + J„+l(x) - J n (x), 

X 



(Bl) 

x 

and Neumann's addition theorem [68] 

oo 

^ J„(xi)J n+k (x 2 ) = J k (xi - x 2 ), (B2) 



where n and k are integer. The addition theorem (B2 1 has two 
special cases: 

oo oo 

2 Jlix) = J] 7„ 2 + ,(x) = 1 (B3) 

fl— — oo n=— oo 

for any integer value of k, and, 

CO 

J] J„(x)J n+k (x) = 0, (B4) 

«=-oo 

for integer # 0. 

The identity J-„(x) = (-l)"J„(x) implies J 2 (x) = J 2 n (x) 
[68 1, and therefore that 



£ n 2m - l J 2 n (x) = 0, 



(B5) 



where j^yi e [0,n — 1], We follow a well-known number 
theoretical approach, which is described in detail in ll66l . The 



number of ways that Eq. (CI i can be satisfied is isomorphic 



to the problem of evaluating the x-independent term in 







2 _n-\\1 
+ X + XT + . . . + 1 



+ X 



n+\\1 



which we multiply by l \ to give 



(1 



r. 

(C2) 



(C3) 



We now require the coefficient of x q(n [) in Eq. (C3 1, which is 
given by 6 1771 



W(2^ 



j(2q\ln(q - j) + q-l\ 



1=0 

'I 



o-M7a 2... 



=Z<-"'( 



_ lv /iV + ,-lV2, 



(C4) 



where we have used (*J = J and set AT = n(g - /). The only 
n-dependent part of W(2q, n) is the binomial coefficient 

IN + q - 1\ _ (N + m - l)(N +m - 2) . . . (N - m + 1) 

\ N-q )~ (2q-l)l ' 

(C5) 

the numerator of which is a polynomial in N of degree 2q - 1 . 
Thus, we may write W(2q,n) = S2 q -\(n) where S2 q -\(n) is a 
polynomial in n of degree 2q - 1 . 



for integer m > 0. In contrast, repeated substitution of the 
recurrence relation ( |B1| | reveals 

oo 2m q oo 

2 n 2m J 2 n (x) = Yj(l) Z J n-<,+2p(x)J n (x), 

n=—oo q=l p=0 n=-oo 

(B6) 

where it is in principle possible, although often tedious, to 



determine the coefficients c qp . With Eqs. (B3l and (B4i we 
can eliminate all terms in Eq. (B6 1 except those where q - 2p. 
Hence, 

oo m 

J] n 2 "'J 2 „(x) = J] a q x 2 " = R 2m (x), (Bl) 

;i=-oo q=\ 

where a q = C2 q , q l2 2q , and we note that the leading-order coef- 
ficient of the polynomial R2m(x) is always a m — y™)/2 2m . In 
particular: 



2n 2 J 2 n (x)= X -, £ ^(*) = ^ + f . (B8) 



APPENDIX C: COUNTING TERMS WHERE G(j f ) = G(j' q ) 

To find the number of terms where G(j q ) = G(j' q ), we re- 
quire that [see Eq. ( |44b| >] 

0'i + h + • • • + h - fi - / 2 -...-/,) = 0, (CD 



APPENDIX D: DERIVATION OF MOMENTS FOR 
GAUSSIAN DISTRIBUTIONS 

1. Second moment 

For a Gaussian initial momentum distribution D(P) = 
exp(-P 2 /2w 2 )/ V27W 2 (w is the standard deviation), the ini- 
tial second-order momentum moment is (P 2 )o = w 2 . Using 



Eqs. (29 1 and ( A22 1, we deduce that, for integer Q = r, 

j2 



(P\ = w 2 + 



2 V2 



l= r 

',7TW 2 J-<* 



dP 



sin (n[a - r]n/2) 
sin 2 ([o' - r]n/2) 



P 2 /2w 2 



(Dl) 

Note that we can replace the f3 in a = ( 1 + 2ji)l with P = k+fi 



without altering Eq. (A22i. We now use Eq. (E3i together 
with 

cos(^[(l + 2P){ - r]n) = {-\f^ cos(_2q€nP), (D2) 



to determine from Eq. (Dl i that 



{ P 2 )„ =w 2 + % + V (-lf^\n - q) 



(D3) 



x f dPcos(2q(nP)e 

*J — oo 



P 2 /2w 2 



' This is equivalent to considering 2q n-sided dice and finding the number of 
ways W(2q,n) of totalling q(n - 1). 
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dPcos(2q£nP)e- p2/2w2 = g-V^W, (D4) 



We now substitute |78 1 

i r 

V2ttw 2 J- 
into Eq. ( |D3| l, which gives the final result: 

(P\ = w 2 + f » + fiYjc-Vf^Xn ~ q)e- 2 " 2eW . (D5) 



9=1 



2. Fourth moment 



The initial fourth-order moment for a Gaussian initial dis- 
tribution is (P 4 )o = 3w 4 . Using a similar approach to that for 
the second-order momentum moment, we use Eqs. p0| and 
( jA22| i to deduce that, for Q, = r, 



(P\ =3w 4 + 



r 

2nw 2 



j p sin (n[a-Q.]n/2) ^_ P 2 /2w 2 



8 V27W 2 J-oo sin*([a - Q]tt/2) 



4>l 



V2tw 2 J— \2 / 
sin 2 (»[o- - Q]tt/2) _ P 2 /2w 2 
sin 2 ([of - Q]tt/2) 



(D6) 



Finally, using <P 2 ) = w 2 , Eq. (|D4]>, and the integral 1781 



Vz 



= f 



dPP l cos(2q£nP)e- pll2wl ^w L {\ - \q L t i n i w L ) 



2»2 J„A 



x e 



•2 9 2 fVw 2 



to simplify Eq. ( D7 1, we find that 



(P 4 ) =3w 4 + ^(2« 2 + 1)« + 3w 2 2 n 
8 

+ ^ ^ ( _ 1)9 (f-r) e -2 f /fV,v 2 
^ 9=1 

i o 2n t 

g J - 2n<? 2 - q + y (2« + 1) 



3*3 



2«-2 



<? 3 2 (12« 2 - 1) 2« 



- 2n</ + 
3 H 3 



tf+y (1-40 



+ 2 ^(-l)^(«-^ W - 2 . 



9=1 
»- 1 



(D8) 



(D9) 



+ 6<p 2 d ^(-l)« (< " r) (n - 4)w 2 (l - 4wV^V) 

9=1 

-2 9 2 /5?r 2 w 2 



x e 



Using Eqs. ( D2 1, (|E3|> and ( E4 1 we deduce from Eq. ( D6 1 that 



APPENDIX E: COSINE EXPANSIONS 



<P 4 > =3w 4 + ^ (2» 2 + 1)« + % + 3<J>§>#! 



30 4 



8 V2ttw 2 
2n 



£(-1)^ 



9=1 



g — 2nq - q 



+ -j (n + D 



f 

.7 -ex 



o #4 2/7-2 

J yc-i) 9(f_r 



dPcos(2q£nP)e- p /2w ~ 

(I2n 2 - 1) 



8 V2^ ^ 



T ^ 2 . 
2nq H 

3 ^ 3 



-q 



+ |(l-4„ 2 ) 



dPcos(2qinP)e 



-P 2 /2w 2 



P 2 - 1 



2(-l)^-')( B -,) 



dPcos(2qinP)e- p /2w 

OO 



60^ 



2 n-1 



J" 



x dPP 2 cos(2q{nP)e 



-P-/2w 2 



(D7) 



In order to perform many of the integrations in Appendix[D] 
we used finite cosine expansions of powers of sin(ni?)/ sin(#). 
To derive such an expansion for sin 2 (n#)/ sin 2 (i?), we first 
consider the summation 

n-1 n-1 n-1 

n + 2 Yj(n ~ q) cos(2qff) =n + n ^] e aq§ + n V 



9=1 



9=1 



9=1 



n-1 



71-1 



9=1 



9=1 



(El) 



The exponential sums in Eq. (El i are geometric sums, or their 
derivatives, which can be evaluated to give, 



n + 2 /fa - q) cos(2<7#) —n + n 



9=1 



2i& _ l 



- n 



1 d e i26n - e m 

2 && e m - 1 
i 



(E2) 



8 g-BOn _ e -2«# 



2 d® e- 2m - 1 
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Differentiating Eq. (E2 1 and identifying a mutual denominator 
then gives the desired result: 



n + 2 /Jjn - q) cos(2qff) 



{e 2i§ - l) 2 



We note that £/„_i(cos(#)) = sin(n)?)/ sin(#) is a Cheby- 
shev polynomial of the second kind, and in this appendix we 
have considered finite expansions of these in terms of Cheby- 
shev polynomials of the first kind, i.e., r„(cos(#)) = cos(n#). 
These are equivalent to discrete Fourier transforms. 



sin (riff) 
sin 2 (#) ' 



(E3) 



Similarly, rewriting cos(2q{P) and converting between terms 
involving powers of q to derivatives of cos(2q§), we find that 



sin 4 (ni?) n 2 



sin 4 (#) 3 



= -(2n z + 1) 



Z 

9=1 
2n-2 

z 

(/=n 



2n 



q 3 -2nq 2 -q + — (2n 2 + 1) 



cos(2gi?) 



1 



2n 



y - 2n<? 2 + (4» 2 - ^)q + ^(1 - 4» 2 ) 



x cos(2g#). 
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